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Abstract We discuss some recent results for the properties of doped antiferromagnets, 
obtained within the planar t-J model mainly by the finite-temperature Lanc- 
zos method, with the emphasis on the comparison with experimental results in 
cuprates. Among the thermodynamic properties the chemical potential and en- 
tropy are considered, as well as their relation to the thermoelectric power. At 
the intermediate doping model results for the optical conductivity, the dynam- 
ical spin structure factor and spectral functions reveal a marginal Fermi-liquid 
behaviour, close to experimental findings. It is shown that the universal form of 
the optical conductivity follows quite generally from the overdamped character 
of single-particle excitations. 



1. INTRODUCTION 

It is by now quite evident through numerous experiments on electronic 
properties that cuprates, being superconductors at high temperatures, are also 
strange metals in the normal phase [1]. On the other hand it also appears 
that most features can be well represented by prototype single-band models of 
correlated electrons, as the Hubbard model and the t-J model. In spite of their 
apparent simplicity these models are notoriously difficult to treat analytically, 
in particular in the most interesting regime of strong correlations. This has led 
to intensive efforts towards numerical approaches [2], mostly using quantum 
Monte Carlo (QMC) and the exact diagonalization (ED) methods. 

The subject of this contribution is the planar t-J model (Hubbard model is 
expected to show similar behaviour in the strong correlation regime), which 
represents layered cuprates as doped antiferromagnets (AFM) and within a 
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single band both mobile charges and spin degrees of freedom, 

H = -tY^ (4^, + H.c.) + jX^(Si • S,- - ^mrij). (1.1) 

Strong correlations are here imposed by strictly forbidding doubly occupied 
sites. So far most calculations were performed for the ground state at T = 0, 
where the standard Lanczos algorithm offers an efficient exact-diagonalization 
analysis of small systems [2]. More recently a novel numerical method, finite- 
temperature Lanczos method (FTLM) [3, 4], has been introduced combining 
the Lanczos method with a random sampling, which allows for an analogous 
treatment of static and dynamic properties at T > 0. 

One of most important conclusion of experimental efforts in the last decade 
is the realization that the electronic phase diagram in the paramater space of 
planar hole concentration Ch and temperature T is quite universal. Materials are 
usually classified as underdoped, optimally doped and overdoped, with respect 
to the highest in a given class. In our analysis we cannot establish the 
superconductivity, so we will use the highest entropy (at low T) as a criterion 
for the optimum doping. In fact both criteria are quite close for real cuprates 
[5] and one could conjecture that this relation is not accidental. In particular 
since in thermodynamic quantities at the same doping also the pseudogap scale 
disappears. 

In this contribution we mainly discuss two topics related to the normal-state 

properties of cuprates. In Sec. II we deal with the thermodynamic quantities: 
entropy, chemical potential and closely related thermoelectric power, all of them 
in the full range of Ch- In Sec.III we concentrate on the appearance and the 
relation between different manifestations of the marginal-Fermi-liquid (MFL) 
behaviour, observed in the optical conductivity, dynamic spin susceptibihty and 
spectral functions at the intermediate (optimum) doping. 

2. THERMODYNAMICS 

Let us first consider thermodynamic quantities, which can be directly derived 
from the grand-canonical sum: free energy density T, chemical potential ji 
and entropy density s. For these the FTLM is particularly simple to implement 
[6, 4], since it requires only a minor generaUzation of the usual Lanczos method. 
Results presented below were obtained mostly for = 20 sites and parameters 
J/t = 0.3, corresponding to the situation in cuprates (where t ~ 0.4eV). 
Note that in small systems only results above certain (size dependent) T are 
meaningful, i.e. in cases below typically T > Tjs ~ 0.1 1. 

Let us first discuss results for s, shown in Fig.(l.l) for various T. It seems 
quite generic feature of such a system that s{ch) exhibits a (rather broad) 
maximum at the intermediate doping ~ c^. The increase of s on doping 
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Figure 1.1 s vs. Ch at several T [4]. For comparison also experimental results for LSCO [5] 
at highest T = 320 /sT ~ 0.07 t are shown. 

can be plausibly related to the frustration between the AFM exchange oc J and 
the hole kinetic energy oc c/jt preferring the FM configuration. This naturally 
leads to a most frustrated situation at ~ J/t. It is quite fortunate that the 
FTLM works best, i.e. Tjg is lowest, just in the cases with large s and large 
frustration, while other methods have difficulties in such situation. E.g. QMC 
is plagued with the minus-sign problem which seems to be intimately related 
to fermionic frustration. 

Even more surprising fact is the magnitude of s at T < J. E.g. at T = 0.1 1 
at we get 40% of s(T = oo). Clearly we are dealing with a system which 
has very low degeneracy (Fermi) temperature Tdeg < J, far below the free 
fermion value T^g^ ~ 8t. Such a conclusion is in agreement with experiments 
in cuprates. In recent years s has been measured in YBCO and LSCO (also 
presented in Fig. (1.1) in a wide doping regime [5] and our results show good 
quantitative agreement. 

For /x/i(r), presented in Fig.(1.2), we mostly do not find a dependence of 
fih at low T, as expected for a normal Fermi liquid, except within the extremely 
overdoped regime c/i > 0.3. In particular, in a broad range 0.05 < Ch < 0.3 we 
find a roughly linear variation //ft (T) = HhiT = 0)-|-a;A;Br, whereby the slope 
a changes the sign at Ch = '-^ 0.15. The variation fj,{ch) at low T has been 
recently deduced experimentally from the shift of photoemission spectroscopy 
spectra in LSCO [7], and the agreement with our results is quite satisfactory 
[4]. From photoemission results as well as from our Fig.(1.2) it is also evident 
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Figure 1.2 Hole chemical potential fXh vs. T at several dopings Ch [4]. 



that ^(c/j < c^) is very flat which would indicate that at T ^ and Ch < 
the compressibiUty k oc —dch/dji is very large or even diverging, as would 
e.g. follow from the phase separation scenario [8] or the singular Cft — > Umit 
[1]. It should be however stressed that the distinction between these scenarios 
could be relevant only at very low T -C J, since both experiment and numerics 
indicate on quite large s, i.e. a distribution over a wide spectrum of states, even 
at T - J/10. 

It is quite helpful to realize that the free energy density J^{ch, T) relates the 
variation of s = dT/dT and ji = dT/dch, i.e. 



ds 
dch 



dT 



dchdT 



(1.2) 



This coimects the maximum s(c^) with the change in slope dnhicD / dT = 0. 
Moreover Eq.(1.2) allows us to discuss more confidently the slope dfi^/dT = 
aks for which we find in the underdoped regime a ~ 2. Although the latter 
has not been so far verified directly for cuprates, one can extract in the same 
regime from the measured ds/dch for LSCO and YBCO at T > 100 K similar 
values a > 1 [5]. It is quite evident that at a > 1 we are not dealing with 
a degenerate Fermi liquid but rather with the nondegenerate doped carriers, 
which is a situation typical for a doped (nondegenerate) semiconductor. One 
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should just recall the standard expression for /x/j in p-type semiconductor, 

c, = P^e'^^^-^^^y^BT ^ = kB\n^ > ks, (1.3) 

where in our notation ~ 1. The constant slope dfih/dT observed in 
our calculations down to T < 0.1 i is a confirmation of such a picture. In 
experimental results for s one should however notice a reduction of a with T, 
but even at T ~ Tc the system is not evidently a normal Fermi Uquid with 
a < 1. 

Another consequence of such a semiconductor picture is an expression for 
the thermopower S, 

b ~ — ~ — ~ abo, (1.4) 

eoT eoT 

where So = kB/eo = 86fiV/K. The validity of this approximation we have 
verified within the t-J model also directly by evaluating the mixed current- 
energy current correlation function and observing that they are proportional 
to the current-current correlation Cjj^j{u) = fj,h{0)C (uj) . The result is in 
Fig.(1.3) is good agreement with the general experimental observation in 
cuprates [9] of a large and rather T-independent S at low doping. In fact 
instead of the usual semiconductor expression for a at c/j -C 1, Eq.(1.4), in a 
strongly correlated system it is more appropriate to use the proper statistics for 
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the t-J model leading to a ln[2(l — c/i)/c;i], which even predicts the change 
of sign at Ch ^ 0.3. 

3. DYNAMICS AT OPTIMUM DOPING 

It has been quite early estabUshed from experiments that cuprates in the 
normal state do not follow the behaviour consistent with the normal Fermi 
liquid. In contrast several static and dynamic quantities at optimum doping can 
be quite well accounted for within the marginal Fermi liquid (MFL) concept 
[10]. Most evident example is the dynamic conductivity ct{u}) which does not 
obey the usual Drude form with a constant rate 1/r but can be well fitted in a 
broad range of oj, T with the generalized MFL form T~^(a;, T) = X{\uj\ + r]T), 
describing also the well established hnear resistivity law p (xT. It has been nat- 
ural to postulate an analogous MFL behaviour for quasiparticle (QP) relaxation 
in spectral functions as e.g. measured by the angle-resolved photoemission 
spectroscopy (ARPES). Only recently, however, the high resolution ARPES 
experiments on BSCCO [11] seem to be in position to confirm beyond doubt 
this behaviour, obeyed in the optimum-doped materials surprisingly even at 
T < Tc for QP along the nodal direction in the Brillouin zone. Most evident 
indication that also spin dynamics follows the MFL concept is the observed 
anomalous NMR and NQR spin-lattice relaxation rate l/Ti{T) ~ const. [12] 
instead of usual Korringa law in metals. 

By calculating using FTLM several related quantities, descibing charge and 
spin dynamics within the t-J model, we established that the MFL concept 
applies well in a broad range of intermediate hole doping 0.1 < c/j < 0.3. 
We discuss here in particular the dynamical conductivity a{u:), the local spin 
susceptibility xl{^) [13] and the QP relaxation rate as obtained from the 
analysis of spectral functions A(k,oj) [14]. Moreover, a^uo) has been been 
found close to a universal form [13], 

1 _ f,~<^/kBT 

a{u) = Co , (1.5) 

UJ 

in a remarkably broad frequency regime < < c<j* ~ 2t, while Co is 
essentially T-independent for T < J. Resulting 17(0; < oj*) is clearly governed 
by T only. Evidently, Eq.(1.5) reproduces the linear resistivity law p = T/Cq 
and is consistent with the MFL scenario for r~^(a;, T), however in a very 
restrictive way since both MFL parameters are essentially fixed. A reasonable 
overall fit can be e.g. achieved by A ~ 0.6 and r] ~ 2.7. When optical 
experiments on (t{u) in cuprates are analysed within the MFL framework quite 
close values for A, rj are in fact reported [15, 16]. In addition, the model results 
reproduce well also the absolute value of and p{T) [4]. 

Analogous universality has been found also in the spin dynamics, in par- 
ticular when looking at the local spin susceptibility XL(a;) and related spin 
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correlation function Sl{uj), 

Xi(a;) = ^tanh(^)5i(a;), (1.6) 



where ^'^(ti;) = ^'^(ti;) + Sl{—oo) is the symmetrized function, having a fixed 
sum rule 

fOO _ 1 

S{u)dw = {{Stf) = -{l-cn). (1.7) 

The most important message of numerical results on spin dynamics at inter- 
mediate doping is that Sl{oj) is nearly T-independent in a broad range of T, 
in particular for T < J. Moreover «S'l(<^) is only weakly doping dependent 
consistent with the sum rule. So we have a conclusion that at intermediate 
doping the more fundamental and universal quantity is the correlation function 
Siiio) and not the susceptibility xl{^), which is the analogy to the relation 
between C{uj) and cr{u!) in Eq.(1.5). For the spin dynamics this is also the 
message of the anomalous NMR 1/Ti in cuprates [12]. As a result x'Li'^)' 
Eq.(1.6), follows the MFL behaviour i.e. at a; < T one gets anomalous T 
dependence x'Li^) ^^1^ ■ 

A MFL-type QP relaxation is extracted within the t-J model also from the 
analysis of the spectral functions ^(k, near the optimum doping [14, 4]. 
For the characterization of QP properties the self energy S(k, cj) is crucial. 
On the other hand the same information can be also expressed in terms of QP 
parameters Zk, Fk, ek. Both definitions are related as 

A(k,a;) = --Im = (1.8) 

^ ' ^ TT a;-S(k,a;) tt (a; - ek)^ + F^ 

For QP near the Fermi surface the hole-part self energy w < is found to 
be of the MFL form, i.e. ImS ~ -7(0; + ^T) with 7 ~ 1.4 and ^ ~ 3.5. 

7 > 1 means an overdamped character of QP, since the the full width at half 
maximum A ~ 2r(e) > e is larger than the QP (binding) energy e. This should 
be contrasted with the electron-like regime a; > where the damping is found 
to be essentially smaller and consequently QP can be underdamped. 

Here we comment on the relation of our results to recent ARPES results in 
BSCCO. The analysis for hole-like excitations in the nodal direction (0, 0) — 
(tt, tt) shows the MFL form with the QP damping F ~ 0.75a; for a; > T and 
r ~ 2.5T for cj < T [11]. This again means an overdamped character of 
hole excitations, since 2r(e) > e. In making the comparison one should take 
into account that F = Z|ImS|. Since at the peak position we find Z ~ 0.5 
experimental and model values appear reasonably close. 

Let us finally discuss the relation of a(yi) and the associated relaxation 
rate 1/r to the QP damping F [17]. In the case of weak scattering one finds 
1/r ~ 2F. In cuprates as well in the t-J model we are apparently dealing with 
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overdamped QP, so the relation is at least questionable. Also, the conductivity 
form Eq.(1.5) appears to be universal, while the QP damping does not seem to 
be parameter free. 

One approach is to approximate the current-current correlation function 
C{u)), which in general replaces Co in Eq.(1.5), by a decoupling in terms of 
spectral functions A(k, u) neglecting possible vertex corrections, i.e. 

Cioo) = ^^K?j du'f{-u')fiu' -u)A{k,uj')A{k,u;' -oj). (1.9) 

In order to reproduce the MFL form of (t{u)) one has to assume the MFL 
form for the spectral function, Eq.(1.8), i.e. the QP damping of the form 
r = 7(|i^| + ^r). We neglect also the k dependence of T and Z. In fact it is 
enough to assume that Tk{uj) is independent of deviations Akj^ perpendicular 
to the Fermi surface. The latter is just what is observed in recent ARPES studies 
of BSCCO [11]. Replacing in Eq.(1.9) the k summation with an integral over 
e with a slowly varying density of states we can derive 

/ deA{e,u')A{e,u;'-u) = — d-lO) 

where T{lo,u;') = T{u}') + T{uj' — lo). We are thus deahng with a function 
C{io) depending only on the ratio u/T, and on MFL parameters 7,^. For 
7 «C 1 we recover via such an analysis C{(jj) strongly peaked at w = and 
consequently MFL-type a{uj) with 1/t{oo) = 2V{uj/2) [10]. No such simple 
relation is valid when one approaches the regime of overdamped QP excitations 
7 ~ 1 or more appropriate 7^ ~ 1. In Fig.(1.4) we show results for several 7 
fixing ^ = TT. For 7 < 0.2 still a pronounced peak shows up at a; ~ 0, on the 
other hand C{ijj) becomes for 7 > 0.3 nearly constant or very slowly varying 
in a broad range of lo/T. 

The main message of the above simple analysis is that for systems with 
overdamped QP excitations the universal form (1.5) describes quite well (t{u) 
for a wide range of parameters. It should be stressed that nearly constant 
C{u < LO*) also means that the current relaxation rate 1/r* is very large, 
l/r* ~ 00* S> 1/r, i.e. much larger than the conductivity relaxation scale 
apparent from Eq.(1.5) where 1/r oc T follows solely from thermodynamics. 

One should also be aware of the upper cutoff scale lu* for the validity of 
the MFL-like QP damping. In the problem considered here it appears that the 
cutoff is directly related to the current relaxation rate u* ~ 1/r* found in the 
t-J model to be extremely high at the intermediate doping, i.e. to* ~ 2t. The 
latter allows for an effective mean free path /* of only few cells, essentially 
independent of T. Such a short I* can be plausibly explained by assuming that 
charge carriers - holes entirely loose the phase coherence in colhsions with each 
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Figure 1.4 Current-current correlation spectra C(a;) vs. uj/T for various 7 at fixed ^ = tt. 



Other due to the randomizing effect of an incoherent spin background. Note 
again that the short correlation length (even at T < Tc) appears also from the 
analysis of ARPES spectral functions A{k,u]) varying Aky along the Fermi 
surface [11]. 

4. DISCUSSION 

Cuprates in their metalUc phase are anomalous in several respects. One im- 
portant conclusion at least for theoreticians is that most of anomalous properties 
are quite well reproduced also in the prototype t-.J model. The analysis of this 
model has been so far restricted to numerical calculations of small systems, 
nevertheless in the T > Tfg window where macroscopic relevance of FTLM 
results is expected the agreement with experiments is even quantitative, without 
any adjustable parameters. Since the behaviour found experimentally is quite 
generic and universal down to lowest T ~ Tc there is no reason to doubt in the 
generality of model results. 

Nevertheless there are open questions of the existence and the origin of low 
energy scales in cuprates as well as in the t-J model. In the underdoped or 
weakly doped regime FTLM shows the indication for the pseudogap scale T*, 
in particular in the uniform susceptibility xo and in the density of states [4]. 
This scale seems to be related to the onset of short range AFM correlations, 
hence T* oc J. Still for T < T* the entropy remains large as manifested by 
experiments and our results. The electron liquid is thus closer to a nondegen- 
erate system of composite particles than to a degenerate Fermi gas. Only at 
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T ^ the entropy is low enough to make the discussion of possible orderings 
or instabilities relevant. 

The origin of the MFL behaviour of several dynamic quantities and of 
the universal form of (t{oj) and XLi^j^) in the intermediate doping has to be 
intimately related to the large degeneracy in this regime. It has been shown 
that QP are essentially overdamped down to lowest T ~ T^. This can be only 
explained by the scattering on spin fluctuations, which mainly contribute to the 
entropy. On the other hand spins are just strongly perturbed by holes introduced 
by doping, so a self-consistent enhancement seems to be the mechanism which 
dominates the relevant physics. Only at low T ~ Tc apparently this behaviour 
breaks down by an emergence of coherence and new energy scales. 
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